Al'tshuler-Aronov correction to the conductivity of 
a large metallic square network 
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We consider the correction Acr ee due to electron-electron interaction to the conductivity of a 
weakly disordered metal (Al'tshuler-Aronov correction). The correction is related to the spectral 
determinant of the Laplace operator. The case of a large square metallic network is considered. The 
variation of Ao ee (LT) as a function of the thermal length Lt is found very similar to the variation 
of the weak localization Actwl (-£<¥>) as a function of the phase coherence length. Our result for 
A(j ec interpolates between the known Id and 2d results, but the interaction parameter entering 
the expression of A<r oe keeps a Id behaviour. Quite surprisingly, the result is very close to the 2d 
logarithmic behaviour already for Lt ~ a/2, where a is the lattice parameter. 

PACS numbers: 73.23.-b ; 73.20.Fz ; 72.15.Rn 



I. INTRODUCTION 

At low temperature, the classical (Drudc) con- 
ductivity of a weakly disordered metal is affected 
by two kinds of quantum corrections : the first one 
is the weak localization (WL) correction, a phase 
coherent contribution that originates from quan- 
tum interferences between reversed electronic tra- 
jectories. This contribution to the averaged con- 
ductivity depends on the phase coherence length 
L v and the magnetic field : Actwl($, Lip)- The 
temperature manifests itself through L v , since 
phase breaking may depend on temperature, e.g. 
if it originates from electron-electron^ or electron- 
phonon 2 ^ interaction. 

In a metal, an electron is not only elastically 
scattered on the disordered potential, but, due to 
the electron-electron interaction, is also scattered 
by the electrostatic potential created by the other 
electrons. At low temperatures, when the elas- 
tic scattering rate (1/T e ) dominates the electron- 
electron scattering rate (1/t cc (T)), the motion of 
the electron is diffusive between scattering events 
with other electrons. In this regime, electron- 
electron interaction is responsible for a small deple- 
tion of the density of states at Fermi energy (called 
the DoS anomaly or the Coulomb dip) and a cor- 
rection to the averaged conductivity as well, the so- 
called Al'tshuler-Aronov (AA) correctional 5 -! 6 -! 7 -^ 
(see RefSi 10 ' 11 ! 12 for a recent discussion). A A and 
WL corrections are of the same order (but this 
latter vanishes in a magnetic field). However, con- 
trary to the WL, the AA correction is not sensitive 
to phase coherence and involves another important 
length scale : the thermal length Lt = \JD/T 
(h = ks = 1). The AA correction, denoted be- 
low Acr 00 (ir), has been measured in metallic wires 
in several experiment er 14 ' 15 ' 16 ! 17 . From the exper- 



imental point of view, AA correction allows to 
study interaction effects in weakly disordered met- 
als, but also furnishes a local probe of temperature 
in order to control Joule heating effect a 15 ' 17 , which 
is crucial in a phase coherent experiment. 

All the works aforementioned refer to the quasi- 
one-dimcnsional (wire) or two-dimensional (plane) 
situations. Quantum transport has also been stud- 
ied in more complex geometries like networks of 
quasi- Id wires. For example several studies of WL 
have been provided on large regular networks in 
honeycomb and square metallic network a 18 ' 19 , in 
square networks etched in a 2DEG2S, and in square 
and dice silver networks^ 1 -. Theoretical studies of 
WL on networks have been initiated by the works 
of Doucot & Ramma l 22 ' 23 and improved by Pas- 
caud & Montambaus^ 4 - who introduced a powerful 
tool^ 5 . : the spectral determinant of the Laplace 
operator, that will be used in the following (see 
also Ref^ 6 -). 

The aim of this paper is to study how the AA 
correction can be computed in networks. In a first 
part we briefly recall how the spectral determinant 
can be used to compute the WL. Then in a second 
part we will consider the AA correction. 



II. SPECTRAL DETERMINANT AND 
WEAK LOCALIZATION 

Interferences of reversed electronic trajec- 
tories are encoded in the Coopcron, solu- 
tion of a diffusion-like equation (dt — D[V — 
2ieA{x)] 2 )V c {x,x';t) = S(x - x')6(t), where A(x) 
is the vector potential. On large regular networks, 
when it is justified to integrate uniformly the 
Cooperon over the network (see Ref^ 7 - for a discus- 
sion of this point) it is meaningful to introduce the 
space-averaged Cooperon V c (t) = J ^V c (x,x;t) 
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where t v = L^/D is the phase coherence time. 
The factor 2 stands for spin degeneracy. We 
have omitted in (|1I2|) a factor l/s where s is the 
cross-section of the wires. The parameter 7 is re- 
lated to the phase coherence length 7 = l/£ 2 
(note that description of the decoherence due to 
electron-electron interaction in networks requires 
a more refined discussio n 28 i 29 ). The spectral de- 
terminant of the Laplace operator is formally de- 
fined as S(j) = dct(7 — A) = Iln(7 + ^n) where 
{E n } is the spectrum of —A [in the presence of 
a magnetic field, A — > (V — 2ieA) 2 ]. The inter- 
est in introducing S(y) is that it can be related 
to the determinant of a V x V-matrix, where V is 
the number of vertices, that encodes all informa- 
tion on the network (topology, length of the wires, 
magnetic field, connection to reservoirs). We la- 
bel vertices by greek letters. l a p designates the 
length of the wire (af3) and Q a p the circulation of 
the vector potential along the wire. The topology 
is encoded in the adjacency matrix : a a p — 1 if 
a and (3 are linked by a wire, a a p = otherwise. 
X a = 00 if a is connected to a reservoir and X a = 
if not. We introduce the matrix 

M a p = 5 a0 ^X a + ^ ^ a<xn coth s/yl ai ^J 

-a a ps/l . . „ (3) 
smh y/jl a/ 3 

where the a a ^ constraints the sum to run over 
neighbouring vertices. Thei 
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S(j)=I[ ^^ detM (4) 

where the product runs over all wires. We now con- 
sider a large square network of size N x x N y made 
of wires of length l a p = a V(a/3). For simplicity we 
impose periodic boundary conditions (topology of 
a torus), which is inessential as soon as the to- 
tal size of the network remains small compared 
to Ltp. At zero magnetic field the spectrum of 
the adjacency matrix is e„ iTO 
2 cos(2mir / N y ) , with n = 1, 
1, 



2cos(2nn/N x ) + 
■ , N x and m = 



, N y . Therefore 
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The calculation of In S (7) involves a sum that can 
be rep 
Using 



be replaced by an integral when N x , N y 3> L v /a 



dxdy 
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K(l/A), (6) 



where K(x) is the complete elliptic integral of first 
kind22, yields^ 
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where the volume of the network is Vol = 2N x N y a. 
We recover the expression of the WL first derived 
by Doucot & Rammal 2 ^. Figure Q] displays the 
dependence of the WL correction as a function of 
the phase coherence length L v . We now discuss 
two limiting cases. 




FIG. 1: Actwl in unit of2e 2 /h as a function of L v /a 
(at zero magnetic field). The dashed line is the Id re- 
sult. The dotted line is the 2d limit eq. fifty . 



Id limit.- In the limit L v <C a (i.e. y^fa ^> 1) 
2e 2 



Act 



WL 
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+ O (e-' 2a / L ^ (8) 



We compare with the result for a wire of length a 
connected at its extremities : Act^ — — tt C^v — 

L 2 

— ). As we can see the dominant terms coincide. 
Deviations appear when L v /a increases since tra- 
jectories begin to feel the topology of the network. 
This is already visible by comparing the second 
terms of the expansions. 

2d limit.- In the limit L v ~^> a (i.e. y^fa <C 1), we 
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Vol^ ln * (7) = ^ 



ln(4L v /a) 



The conductivity reads 
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with C W l = ^ + 1 - 0.608 . As noticed in the 
beginning of the section, eqs. (|8|10|) should be di- 
vided by the cross-section s of the wires. In the 2d 
limit, diffusive trajectories expand over distances 
larger than a and feel the two dimensional nature 
of the system, being the reason why (|10|) is reminis- 
cent of the 2d result. It is interesting to point that 
the network provides a natural cutoff (the length 
of the wires, a) while the computation of the WL 
for a plane in the diffusion approximation requires 
to introduce a cutoff by hand for lower times in 
eq. ([1]), which is the elastic scattering time r e . In 
this latter case the constant added to the logarith- 
mic behaviour is not well controlled since it de- 
pends on the cutoff procedure (the computation of 
the constant for a plane requires to go beyond the 
diffusion approximation and leads to^L Act£ ane — 



since 
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AL'TSHULER-ARONOV 
CORRECTION 



At first order in the electron-electron interac- 
tion, the exchange term is the dominant contribu- 
tion to the correction to the conductivit y 8 ' 1 ' 1 1 ' 3 2 
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where U(q,u>) is the dynamically screened inter- 
action. Within the RPA approximation and in 
the small q and lu limit, the interaction takes the 
forrn^ 3 - U (q, lu) ~ " u ^? <? where po is the den- 
sity of states per spin channel. Replacing the 
Drude conductivity by its expression ctq = 2e 2 poD 
and performing an integration by parts, we get 
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After Fourier transform, the result can be cast in 
the formii : 
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(13) 

For the exchange term considered here, one finds 
X a = 4/d. Further calculation yields 8 - \ a ~ 
| — |F, where F is the average of the interaction 
on the Fermi surface (sec definition in Refs4£). 
This expression of A CT is valid in the pcrturba- 
tive regime, F <C 1 ; nonperturbative expression 
is given in Refs.— Vdit) is the space inte- 
grated return probability Vd{t) = / ^Vd(x,x;t), 
where Vd{x, x 1 ; t) is solution of a classical diffusion 
equation similar to the equation for V c (x,x';t), 
apart that it does not feel the magnetic field : 
[d t -DA]Vd(x,x';t) = S(x-x')S(t). Therefore the 
Laplace transform of Vd(t) is given by 9 7 ln5(7) 
with 6 a j3 = 0. It is interesting to point out that 
([13")) has a similar structure to (P} with a different 
cutoff procedure for large time. It also involves 
a different scale : the temperature dependence of 
Acr ee is driven by the length scale Lt instead of 
L v for the weak-localization correction A<7\vl- 



Up to eq. (|13p the discussion is rather general 
and nothing has been specified on the system. 
We have seen in section [XT] that the WL for the 
square network presents a dimensional crossover 
from Id to 2d by tuning L v /a. A similar dimen- 
sional crossover occurs for the AA correction by 
tuning Lt/cl as we will see. This remark raises the 
question of the dimension d in eq. (fTTj) . To answer 
this question we should return to the origin of the 
factor l/d : the current lines in the conductivity 
(Xij produce a factor q^q-j replaced by &ij\q 2 after 
angular integration. Since in a network the diffu- 
sion in the wires has a Id structure (provided that 
W <C L T ~ ^D/lu, where W is the width of the 
wires), the dimension in A CT is d = 1. Therefore we 
have for the network A£ ctwork ~ 4 - \F. 

If one now expands the thermal function in (|13[) 
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we can also relate Aa ee to the spectral determi- 
nant. We obtain : 



oc 1 
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(15) 

which is the central result of this paper. It is the 
starting point of the discussion below. 

Application to the case of the square network. - We 
have to compute 7 2 ^p- IhS^t). We start from ([7]) 
and compute its second derivative. We obtain after 
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some algebra : 

2 00 1 

a \ e a sr^ 1 

Acr cc = -A CT — - > - 1 
h 8 ^— ' " 



where the function ip(x) is given by 
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E(x) being the complete elliptic integral of second 
kind». The function ip(x) is plotted in figure [2] 
and its limiting behaviours are easily obtained 3 ^ : 

(p(x) = -+0(x 2 ) fora-^0 (18) 

= - - 4 + 0(xe- 2x ) for x 00 (19) 
a; 

The Lt dependence of AA correction on a square 
network is displayed on figure [31 where we have 
plotted Aa ee (LT) given by eq. fl"6]) . The di- 
mensional crossover now occurs by tuning the ra- 
tio Lt/ci. We consider the two limits. 
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FIG. 3: The continuous line IS l\0 'cc in unit of Ao- 77 
as a function of Lt/cl (the series H16[) is computed nu- 
merically). The dashed line is the Id limit, eq. \2U[) . 
and the dotted curve is the 2d limit, eq. $21\) . 



haviours of f(x) that the first sum diverges loga- 
rithmically with AT while the second brings a neg- 
ligible contribution of order Af° . Therefore : 



FIG. 2: The function ip(x) of eq. (H\) . 

Id limit.- For Lt <C a we can replace the expan- 
sion (1191) in the scries ([TBI). Therefore 
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with 3 -^2X 



0.782. The dominant term again 

coincides with the one for a connected wire ^ 10 ' 11 
while the second differs by a factor 2, as for the WL 
[see discussion after eq. (jHJ]. 

2d limit.- In the limit Lt 3> a we introduce 
Af = (Lt/cl) 2 and cut the sum (fTBj) in two pieces : 
J2T = Ei^ + Ea'- ^ i s clear from the limit be- 



A(7 n 



-Act — a 
h 



- ln(£ T /a) + C" , 

7T 



(21) 



The constant is estimated numerically. We find 
C cc ^ 0.56. 



The two eqs. (|20|21j) should be divided by the 
cross-section s of the wires. 

The two functions Actwl^ = 0, L v ) (figure Q} 
and Aa ee (LT) (figure [3|) arc very similar. Apart 
from the prefactors 2e 2 /h and X a e 2 /h which ac- 
count respectively for the spin degeneracy and the 
interaction strength, the linear behaviours at the 
origin have a different slope (1 and 0.782) and 
the logarithmic behaviours are slightly shifted : 
Cwl ^ 0.61 and C cc ~ 0.56. 
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IV. COMPARISON WITH EXPERIMENTS 

The AA correction has been recently measured 
by Mallet et ai^i in networks of silver wires with 
3 10 4 and 10 5 cells, lattice spacing a = 0.64 fim 
and diffusion constant D ~ 100cm 2 /s. The dif- 
fusion constant D has been measured separately 
(through measurement of the Drude conductance) , 
therefore we can compare our result (|16p with ex- 
periment using one fitting parameter only : the 
interaction parameter X a . The 2d logarithmic 
behaviour (f2Tj) has been observed in the range 
100 mK< T < 1 K from which the value A° xp ~ 3.1 
was extracted, in agreement with similar mea- 
surements performed on a long silver wire for 
whichSi^ A« x p< wiro ~ 3.2. We now compare with 
the theoretical value : for silver Fermi wavelength 
is kp 1 = 0.083 nm and Thomas-Fermi screening 
length = 1/ \J%-KpQe? = 0.055 nm. In the 
Thomas-Fermi approximation, the parameter F is 
given byii F = (5^) 2 ln[l + (^f-) 2 }, therefore 
F ~ 0.58. Using the Id nonperturbative expres- 
sion^ X a = 4 + §(y/l + F/2 - 1 - F/4), we get 
A™ 1 ~ 3.24, close to the experimental value. 

V. CONCLUSION 

Equations (jl5|16| are our main results. The first 
one shows that AA and WL can be formally re- 
lated : 

Aa cc (L T ) = ^ (22) 

00 1 r d 2 

L f> T 

The validity of this relation is the same as for 
eqs. (|l|2p : the system should be such that it 
is meaningful to average uniformly the nonlocal 
conductivity a{r,r') to get the local conductivity 

a = J d vo[ a ( r > r ')- A similar discussion has been 
proposed to relate WL and conductivity fluctua- 
tions (see appendix E of Rcf.— ). 

Our starting point (fTTj) is a formulation in the 
Fourier space, what implicitly assumes translation 
invariance. Whereas this assumption seems rea- 
sonable for a large regular network such as the 



square network studied in this article, its valid- 
ity is not clear for networks of arbitrary topology, 
what would need further developments. 

We have computed the AA correction in a large 
square network and shown that the result inter- 
polates between the Id, eq. ([20|) . and a 2d result, 
eq. (|2"Tj) . Interestingly, the 2d limit in a network 
involves a Id constant A£ ctwork ~ 4 - §F, what 
is confirmed by experiments, as discussed in sec- 
tion nvi 

The interest of the network compared to the 
plane is to control the constant C cc of eq. (f2"Tj) : 
for a plane, a cutoff must be introduced in eq. (|13p 
at short time t ~ r e and the constant C ee is 
replaced by a number that depends on the pre- 
cise cutoff procedure. Experimentally, it would 
be interesting to observe the crossover from (f20|) 
to (|2ip by varying Lt/cl. This was not pos- 
sible in experiments of Mallet et a£i described 
in section IIVI because measurements are compli- 
cated by the fact that electron-phonon interaction 
also brings a temperature-dependent contribution, 
A<7 e _ p h, at high temperature (above few Kelvins). 
The conductivity is given by a = ao + Actwl + 
Acr ec + Adc-ph. The WL can be suppressed by a 
magnetic field however the electron-phonon contri- 
bution is difficult to separate from Aer co . There- 
fore the network should be patterned in a way such 
that the crossover ld-2d remains below T ~ 1 K 
where Acr _ p h is negligible. As an example we con- 
sider the silver networks studied in Ref^i for which 
L T = 0.27 x T- 1 / 2 (L T in /im and T in K). In 
order to see clearly the Id and the 2d regimes it 
would be convenient to study two networks with 
different lattice spacings. If temperature is con- 
strained by 10 mK< T < 1 K, for a = 0.5 /im we 
have 0.5 < Lt/cl < 5, which probes the 2d regime 
over one decade. A second lattice with a ~ 5 /im 
would allow to probe the Id regime since in this 
case 0.05 < L T /a < 0.5. 
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